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ABSTRACT 

in ; 

O . We used HST/STIS to obtain the spectrum of molecular hydrogen associated 

^ ! with the damped Lya system at z^hs = 1-7765 toward the quasar Q1331+170 

^ ! at Zem = 2.084. Strong H2 absorption was detected, with a total H2 column 

, ^ ! density of A^(H2) = (4.45 ± 0.36) x 10^^ cm~^. The molecular hydrogen fraction 

! is = N^-^^2Nn ~ (^'^ -'- 0-''')%' which is the greatest value reported so far 

i in any redshifted damped Lya system. This results from the combined effect 

^ i of a relatively high dust-to-gas ratio, a low gas temperature, and an extremely 

^ ! low ambient UV radiation field. Based on the observed population of J states, 

! we estimate the photo-absorption rate to be i?abs = (7.6 ± 2.4) x 10"^^ s~^, 

^ ! corresponding to a local UV radiation field of J(IOOOA) ^ 2.1 x IO^^-^^^qqq^ ^5 

■ where J ° is the UV intensity at 1000 A in the solar neighborhood. This 

10 . 1OOOA,0 ^ 

O ■ is comparable with the metagalactic UV background intensity at this redshift, 

1^ • and implies an extremely low star formation rate in the absorber's environment. 
Qh' 

I ■ We construct a simple model to describe the structure of the H2 absorber, with 

2 ' a best-fit total hydrogen number density of n(H) 0.2 cm~^ and an electron 

^ ■ temperature of ~ 140 K. Assuming spherical symmetry, the mass of the H2 

^ ' cloud is estimated to be ~ 6.5 x 10'' Mq, larger than the masses of most giant 

v~j ■ molecular clouds (GMCs) in the Milky Way and nearby galaxies. The extinction 

■ of Q1331-I-170 due to the intervening DLA is -Eb-v ~ 0.037, and we also find 



that the extinction by DLAs with firm H2 detections is significantly greater than 
those for which only upper limits of fu^ have been put. The observed CO-to-H2 
column density ratio is < 2.5 x 10~^, which is similar to the value measured 
for diffuse molecular clouds in the Galactic ISM. Finally, applying the inferred 
physical conditions to the observed C I fine structure excitation (Songaila et 
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al. 1994), we estimate the cosmic microwave background temperature to be 
TcMB = (7.2±0.8)K at z = 1.77654, consistent with the predicted value of 
7.566 K from the standard cosmology. 

Subject headings: quasars: absorption lines — quasars: individual (Q1331+170) 
— ISM: molecules 



1. Introduction 

Damped Lya absorbers (hereafter DLAs) are defined as quasar absorption systems with 
neutral hydrogen column density, iV(Hl) > 2 x 10^° cm~^. They are generally considered to 
be progenitors of present-day galaxies, in the form of fast rotating disks (Prochaska & Wolfe 
1997), or merging protogalactic clumps (Haehnelt et al. 1998). They serve as an important 
gas reservoir for star formation at high redshifts (e.g. Storrie-Lombardi & Wolfe 2000). 

Molecular gas, especially molecular hydrogen, is an important ingredient of star forma- 
tion, since stars form by cooling via molecular lines. The search for H2 in DLAs can be 
carried out by observing the H2 absorption lines in the Lyman X^E^ — ^B^S+ and Werner 
X^S+ — *^C^E^ bands, with rest wavelengths from ~ 1220A to the Lyman hmit (Abgrall 
et al. 1993a, 1993b). Ge & Bechtold (1999) conducted a survey for H2 absorption in 13 
DLAs, using the MMT Spectrograph. Two have been confirmed with significant H2 fraction; 
upper hmits in the range of 10~^ to 10~^ were put on the H2 fraction in other systems (Ge 
&; Bechtold 1997, 1999, Ge et al. 2001). More recently, based on the high resolution spectra 
obtained with UVES at the ESO VLT, Lcdoux et al. (2003) searched for H2 in a sample of 
33 DLAs, with firm detections of H2 absorption in 8 of them. By contrast, almost every line 
of sight through Galactic or Magellanic gas clouds with comparable H I column density has 
strong H2 absorption (e.g. Shull et al. 2000). 

In the interstellar medium, H2 molecules are expected to form on the surface of dust 
grains by physical or chemical adsorption (Cazaux & Tielens 2002). Therefore one contribut- 
ing factor to the lack of H2 in DLAs is the low dust-to-gas ratios in these systems (e.g. Pei et 
al. 1991, Pettini et al. 1994, Vladilo 1998, Prochaska & Wolfe 2002), perhaps the result of 
the fact that quasar surveys tend to select objects which are bright and blue (e.g. Gregg et 
al. 2002, Khare et al. 2004). A second contributing factor to the lack of H2 may be the high 
gas temperature in DLAs which makes H2 formation onto dust grains inefficient (Petitjean 
et al. 2000, Liszt 2002). This is consistent with the suggestion that unlike the local ISM, 
the gas in DLAs might be predominantly warm (Lu et al. 1996, Prochaska & Wolfe 1999). 
Finally, H2 molecules are easily dissociated by UV photons in the energy range of 11.3 — 
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13.6 eV. Therefore the low H2 fraction in DLAs may be the result of a strong UV radiation 
field (e.g. Black et al. 1987), which could in turn be a signature of high star formation 
rate (Wolfe et al. 2003a, 2003b, 2004). Recently, Hirashita et al. (2003) showed that the 
deficiency of H2 in DLAs may also be a result of the expected low volume filling factor of 
H2 in environments with a strong UV radiation field and low dust-to-gas ratio, rather than 
a real absence of H2. 

Q1331+170 (at z^m = 2.084) is one of the first high redshift quasars discovered (Baldwin 
et al. 1973), and has a damped Lya absorber aX z = 1.7765 (Carswell et al. 1975) with a 
neutral hydrogen column density of A^hi = 1-5 x lO^^cm"^ (Chaffee et al. 1988, Prochaska & 
Wolfe 1999). The 21-cm absorption line was also detected (Wolfe & Davis 1979, Chengalur 
& Kanekar 2000), with a spin temperature of Tg ^ 770 K. This implies the presence of a cold 
neutral medium (CNM) phase conducive to the efficient formation of H2 onto dust grains 
(Petitjean et al. 2000). One of the most interesting features in the absorption spectrum of 
Q1331+170 is the presence of strong C I and C I* multiple lines near A ~ 1656 A (Meyer et 
al. 1986, Songaila et al. 1994). C 1 absorption has been found to be rare in quasar absorbers 
as compared to Galactic diffuse clouds as a result of higher UV radiation fields, or lower dust 
abundances (Ge et al. 1997, 2001). Since these conditions would also be inhospitable to an 
appreciable formation of H2, the damped Lya absorber toward Q1331+170 is an excellent 
candidate for a search for H2 absorption. 

Detecting C 1 is also of cosmological interest, since it is insensitive to the local physical 
conditions, but depends strongly on the cosmic background level (Liszt 2002). Therefore the 
fine structure transitions of C I can be used to measure the CMB temperature as a function 
of redshift (e.g. Bahcall et al. 1973). In the case of the DLA toward Q1331-I-170, the 
excitation temperature of neutral carbon determined by Songaila et al. (1994) is 7.4 ± 0.8 K 
for one of two velocity components at z — 1.77654, in agreement with the predicted value 
of TcMB = 7.566 K at that redshift. This result implies that local excitation of C I is 
negligible. Since the population of H2 at different J states can also be used to determine 
physical parameters, such as kinetic temperature, neutral hydrogen number density, and UV 
radiation field, observations of molecular hydrogen and its rotational excitation provide an 
independent verification of the observed C I excitation. 

In this paper, we present our detection of redshifted H2 in the high resolution absorp- 
tion spectrum of Q1331+170 taken with the Space Telescope Imaging Spectrograph (STIS) 
aboard the Hubble Space Telescope (HST). Sec. 2 describes the observations and basic data 
reductions. In Sec. 3, we describe the Voigt profile fitting of H2 lines, and present the best-fit 
physical parameters. Discussions are presented in Sec. 4, and conclusions in Sec. 5. 
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2. Observations and continuum fitting 

In HST GO programs 7271 and 9172, we used HST/STIS to obtain the near UV spec- 
trum of Q1331 + 170 in the wavelength range of 2500A to 3120A. The spectrum was obtained 
with the E230M grating with 0.2 x 0.2 aperture, and the NUV/MAMA detector in the 
ACCUM operating mode, resulting in spectral resolution oi R — 30, 000 (or 10 km s~^) at 
A ~ 2700A. The total exposure time was 52,298 seconds, yielding spectra with S/N~7 above 
the Lyman limit of the damped Lya absorber (~2530A). The data sets are hsted in Table 
1. 

The data were pre-processed with the IRAF cahbration routine CALSTIS (v2.7) for 
nonlinearity correction, dark subtraction, flat-fielding, wavelength calibration, and extraction 
of the 1-D spectrum. Rebinning was performed on the raw images for Nyquist sampling (i.e. 
2 pixels per resolution element). Instead of using automatic wavelength calibration files 

taken for each Mode Selection Mechanism (MSM) setting, we obtained additional exposures 
of comparison lamps during each orbit. Typical observing conditions may cause thermal 
drifts of 0.1 pixel per hour (see Sec. 11.2 of Space Telescope Imaging Spectrograph Instrument 
Handbook for Cycle ISioi more details). For our longest single exposure, this corresponds to 
a wavelength uncertainty of ~ 0.007A at A ~ 2700A introduced by wavelength cahbration, 
or a velocity uncertainty of ~ 0.8 km s~^. 

The entire flux-calibrated spectrum of Q1331-I-170 is shown in Fig. 1, smoothed by 5 
pixels. The continuum was determined by a 4th order polynomial fit using the data longward 
of the Lyman hmit of the DLA (fs 2530A) . Regions with strong absorption features identified 
by eye were excluded. The dotted hue in Fig. 1 shows the 1 a uncertainties of the observed 
spectrum. The continuum fit suffers from considerable uncertainties, due to the numerous 
absorption features in the Lya forest, as well as H2 absorption lines in the Lyman and Werner 
bands. For this reason, we checked our continuum fit according to the following prescription. 
We separate the whole spectrum into several continuous wavelength regions, and for each 
region, we construct the histogram of all the data points with normalized flux greater than 
1. Ideally, the histogram can be well represented by a (half-) Gaussian distribution with 
a standard deviation corresponding to the normalized 1 a error of our spectrum. Next, we 
adjust the best-fit continuum in the region under consideration until the calculated Gaussian 
standard deviation is approximately equal to the average normalized 1 a error in the same 
region. The above procedure was performed in an iterative manner, since an adjustment of 
the continuum also changes the (normalized) 1 a error used for comparison. This procedure 
can be used to determine both the overall amplitude and shape of the continuum. 

We estimate the accuracy of continuum fitting by calculating the reduced x^-deviation 
between the probability distribution function (PDF) of normahzed flux adopting the best-fit 
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Table 1. Log of HST/STIS observations 



Dataset 


Program 


UT date of observation 


Exposure (sec) 


o56m0101 


7271 


1999-08-20 


2076 


o56m0102 


7271 


1999-08-20 


5794 


ODomUiUl 


7271 


1 c\c\c\ no o 1 

1999-08-21 


2076 


o56m0202 


7271 


1999-08-21 


5794 


o6cy0101 


9172 


2002-07-08 


2097 


o6cy0102 


9172 


2002-07-08 


5794 


o6cy0201 


9172 


2002-08-12 


2097 


o6cy0202 


9172 


2002-08-12 


5794 


o6cy0301 


9172 


2002-08-19 


2097 


o6cy0302 


9172 


2002-08-19 


5794 


o6cy0401 


9172 


2002-08-19 


2097 


o6cy0402 


9172 


2002-08-19 


5794 


o6cy0501 


9172 


2003-03-16 


2097 


o6cy0502 


9172 


2003-03-16 


5794 


o6cy0601 


9172 


2003-03-12 


2097 


o6cy0602 


9172 


2003-03-12 


5794 


o6cy0701 


9172 


2002-07-10 


2097 


o6cy0702 


9172 


2002-07-10 


2897 
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Ql 331 + 170 STIS E230M 




2500 2'60fi 2700 2800 2900 3000 
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Fig. 1.— The STIS E230M spectrum of Q1331+170. The best-fit continuum (the thick 
sohd fine) is superposed on the fiux-cahbrated spectrum. The dotted fine corresponds to 
1 (J uncertainties of the observed spectrum. The data have been smoothed by 5 pixels. The 
two dashed fines represent tfie best-fit continuum scaled 15% lower and 15% higher. Also 
shown on the figure are the synthetic profiles of Lyman series lines, including Ly/?, Ly7, Ly5 
and Lye. The synthetic profiles are calculated with a neutral hydrogen column density of 
1.5 X 10^^ cm~^ and a 6-parameter of 19.6 km s~^ (see text). 
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continuum and that with continuum level scaled by a certain factor (ranging from 50% lower 
to 50% higher). Here only pixels with normalized flux greater than 1 are considered. This 
calculation shows that the best-fit continuum is accurate within 15% at 83% confidence level, 
and accurate within 25% at 98% confidence level, based on the distribution. 

The two dashed hues in Fig. 1 correspond to the best-fit continuum scaled 15% lower 
and 15% higher. This represents a very conservative estimate of the uncertainty in the 
continuum fit; clearly the dashed lines are inconsistent with the expected continuum level in 
portions of the spectrum. The effects of the possible uncertainties in continuum fitting on 
our final results are discussed in Sec. 3.3. 



3. Analysis of molecular hydrogen lines 

3.1. Detection of molecular hydrogen 

Previous spectroscopic studies show that there are at least five redshift systems with 
detectable metal absorption lines along the fine of sight toward Q1331-I-170. They include 
the damped Lya system a,t z — 1.7765, as well as absorbers Bit z — 0.7441 (Sargent et al. 
1988), z = 1.3284 (Steidel & Sargent 1992), z = 1.4458 (Sargent et al. 1988) and z = 1.7864 
(Lanzetta et al. 1987). 

We searched for molecular hydrogen absorption associated with the z — 1.7765 DLA and 
two other absorbers z — 1.4458 and z — 1.7864, since their redshifts allow considerable 
overlap between the wavelength coverage of our spectrum and the redshifted Lyman and 
Werner bands of H2. By comparing with the synthetic H2 spectrum, we detected strong 
molecular hydrogen absorption associated with the z = 1.7765 damped Lya system, whereas 
the H2 line pattern does not match the observed spectrum for the other two absorbers. 

For the damped Lya system a,t z — 1.7765, we calculated the significance level of H2 
fines for each J state. We define the significance level, SL associated with state J as 

_ E J[l - /.(A)]ciA 

[Ej^f(A)rfA]V2' U 

where /i(A) and o"j(A) arc the normalized flux and 1 a error within an individual H2 line, and 
the summation is over all the selected H2 lines for the same J state. The results show that 
based on our spectrum, J = 0, 1, 2, 3, 4 fines are detectable at 4.4 cr, 5.1 cr, 9.0 cr, 8.9 a and 
3.7 a significance level, respectively. The J > 5 lines were not detected at 3 cr significance 
level. The number of unblended H2 fines included in the above calculation is two for J — 0, 
three for J — 1, six for J — 2, ten for J = 3, two for J = 4, and five for J — 5, respectively 



-8- 



(see also Table 2). Therefore, we conclude that we detect molecular hydrogen absorption 
from the J — to J — A states, associated with the damped Lya system at 2; = 1.7765. 

In Fig. 2, we mark the expected positions of metal lines from all the known absorbers, 
including the Milky Way interstellar medium, on the observed spectrum (solid line) with 1 a 
errors (dotted line). The atomic line hst was taken from Table 2 of Prochaska et al. (2001). 
Also marked on Fig. 2 arc all the H2 lines from J = to 4 in the Lyman and Werner bands 
for the z — 1.7765 absorber. 



3.2. Voigt profile fitting of molecular hydrogen lines 

We used VPFIT (version 5) to simultaneously fit Voigt profiles for all the H2 hues. 
All lines with J = 0, 1, 2, 3 and 4 were fit simultaneously, assuming one fo-parameter and 
redshift for all lines and different column densities for each J state. These free parameters 
were adjusted independently to minimize the value of the overall fitting. We adopted 
the STIS instrument profile for the echelle E230M and 0.2 x 0.2 aperture. Those H2 lines 
with large reduced deviations {xl > 1-2) were excluded, to avoid possible cases in which 
H2 absorption may be contaminated by Lya forest lines. This procedure was iterated until 
no more lines were rejected. Our final model includes 20 different regions with 26 molecular 
hydrogen fines for simultaneous fitting (two J — fines, three J — 1 fines, seven J — 2 
fines, twelve J — 3 fines and two J — A fines). We also note that the physical parameters 
determined from Voigt profile fitting (see below) are insensitive to our adopted maximum 
reduced deviation of 1.2. 

Table 2 provides the results for all the H2 fines included in the final fitting. Rest-frame 
wavelengths, /-values and damping parameters F are taken from Abgrall et al. (1993a, 
1993b). The sixth column gives the reduced value for all the 20 fitted regions. We show 
the best-fit Voigt profiles for different H2 lines in Fig. 3 (solid lines), superposed on the 
observed spectrum with 1 a errors. The Voigt profiles for all the H2 lines with J = to 4 in 
the Lyman and Werner bands are also shown by the dashed line in Fig. 2. 



Fig. 2a.— The STIS E230M spectrum of Q1331+170 (solid line) with line identifications 
and la errors (dotted line). The molecular hydrogen lines from J = to 4 in the Lyman and 
Werner bands, as well as the expected positions of metal lines from known metal absorbers 
are marked. Numbers in the parenthesis following the metal line identifications indicate 
different absorption systems: - the damped Lya absorber at z = 1.7765, 1 - the absorber 
at z = 0.7441, 2 - the absorber at z = 1.3284, 3 - the absorber at z = 1.4458, 4 - the absorber 
at z = 1.7864, and 5 - the Milky Way. The dashed line gives the H2 synthetic spectrum, 
calculated with the best-fit parameters in Table 3. 
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Fig. 2b. — Same as Fig. 2a. 




Fig. 2c. — Same as Fig. 2a. 




Fig. 2d. — Same as Fig. 2a. 




Fig. 2e. — Same as Fig. 
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Fig. 2f. — Same as Fig. 2a. 




Fig. 2g. — Same as Fig. 
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Fig. 2h. — Same as Fig. 2a. 
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Fig. 3. — The best-fit Voigt profiles for the 26 H2 hnes used in our simultaneous fitting (solid 
lines), overplotted on the observed profiles with 1 a errors. 
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Table 2. H2 absorption lines associated with the DLA toward Q1331+170 



No. 


Identification 




/ (10-^) 


r (109 s-i) 




1 


L17P2 


927.02 


0.233 


0.804 


1.01 


2 


L16P1 


932.27 


0.240 


0.835 


0.32 


3 


W4P3 


934.79 


0.709 


1.101 


0.90 


4 


L16R3 


935.57 


0.657 


0.830 


0.80 


5 


L15R3 


942.96 


0.574 


0.862 


1.06 


6 


L13P3 


960.45 


0.493 


0.939 


0.73 


7 


L10P3 


987.77 


0.893 


1.070 


1.10 




W1Q2 


987.97 


3.650 


1.157 




8 


L9P3 


997.82 


0.876 


1.120 


0.76 


9 


L8P2 


1005.39 


0.991 


1.177 


0.46 


10 


L8R3 


1006.41 


1.572 


1.167 


0.73 


11 


L7P2 


1016.46 


1.023 


1.236 


0.72 


12 


W0P4 


1017.38 


0.995 


1.178 


1.10 




L7R3 


1017.42 


1.838 


1.224 








iUDO.ZO 


n nno 

u.yuz 


i.44o 


U. lO 


14 


L4R3 


1053.98 


1.336 


1.434 


0.85 


15 


L4P3 


1056.47 


0.956 


1.445 


0.59 


16 


L3P4 


1074.31 


0.774 


1.525 


0.98 


17 


L2R0 


1077.14 


1.166 


1.630 


0.62 




L2R1 


1077.70 


0.769 


1.626 




18 


L2P2 


1081.27 


0.470 


1.630 


0.93 




L2R3 


1081.71 


0.636 


1.613 




19 


LIRO 


1092.19 


0.578 


1.738 


0.88 




LlRl 


1092.73 


0.378 


1.734 




20 


L1P2 


1096.44 


0.237 


1.738 


0.94 




L1R3 


1096.73 


0.305 


1.719 
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The best- fit results of H2 column densities, 6-parameter and redshift are listed in Table 3. 
The 1 a errors were derived from ^ fitting by VPFIT, and are mainly dependent on photon 
statistics. The total molecular hydrogen column density is iYnj = (4.45 ± 0.36) x 10^^ cm~^. 
We can place an upper limit on the J = 5 column density as A^j=5 < 6.3 x 10^'^ cm"^, at 2 o" 
significance level. The 1 a error in redshift given in Table 3 corresponds to the uncertainty in 
line centroiding, which is equivalent with 0.9kms~^ in velocity space. We mentioned in Sec. 
2 that the velocity uncertainty caused by wavelength calibration is 0.8 km s~^. Therefore the 
total 1 a error is 1.2 km s~^, which is about 10% of the spectral resolution. 

The molecular hydrogen fraction is defined as 

TVhi + 27Vh, ■ 

Adopting a neutral hydrogen column density of A^hi = (1.50 ± 0.14) x 10^^ cm~^ (Prochaska 
& Wolfe 1999), we derived a molecular hydrogen fraction of /hj = (5.6 ± 0.7)%. This is 
the largest value reported so far in any redshifted damped Lya system. The profiles of the 
Lyman scries lines are shown in Fig. 1, assuming A^hi = 1-50 x 10^^ cm^^ and h = 19.6km s^^. 
The 6-parameter is directly scaled from the 6-parameter determined from Voigt profile fitting 
of molecular hydrogen lines, taking into account the difference in particle mass between the 
two species. The synthetic hue profiles are in broad agreement with the observations. The 
fact that H I lines are heavily blended with other H2 lines makes Voigt profile fitting to 
derive A^hi difficult. Therefore we simply adopt the neutral hydrogen column density from 
previous work (Prochaska & Wolfe 1999). 



3.3. Uncertainties due to curve-of-growth and continuum fitting 

The uncertainties in the measured H2 column densities at different J states depend on 
the location of the corresponding lines on the curve of growth. We show in Fig. 4 the curve 
of growth for the 26 H2 hues used in our simultaneous Voigt profile fitting. Note that the 
equivalent widths arc calculated with the best-fit parameters shown in Table 3, since some 
of the lines are blended and therefore not appropriate for direct measurements of equivalent 
widths from the observed spectrum. Lines with different J values are marked with different 
symbols. The background dots in Fig. 4 bracket the full range of the curve of growth for H2 
lines, taking into account the variations in the atomic data (Abgrall et al. 1993a, 1993b). 
The figure shows that the J = 0, 1 hues he on the damped part of the curve of growth, while 
J = 2, 3 lines on the flat part, and J — A hues on the hnear part of the curve of growth. 
Since the determination of column densities is very sensitive to the choice of 6-parameter 
on the flat part of the curve of growth, the uncertainties in A^j=2 and A^j=3 are relatively 
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larger than the column densities at lower J states. Although J — 4 lines are on the hnear 
part of the curve of growth, the 1 a error in Nj=4 is also relatively large due to the low 
signal-to-noise ratio for these weaker lines. 

Uncertainties in column density are also likely introduced by uncertainties in continuum 
fitting. To investigate this effect, we scaled our continuum by a factor ranging from 25% 
lower to 25% higher than the best-fit continuum. Then we did simultaneous Voigt profile 
fitting in the same way as described in Sec. 3.2. The results are shown in Fig. 5, for the 
measured column densities at different J states. We find that the errors in column density 
at all J states caused by 1 cr uncertainties in continuum fit arc considerably smaller than 
the 1 a error in column density derived from Voigt profile fitting. In the extreme case that 
if our best-fit continuum is under-estimated by a factor of 25% (which can be excluded by 
98% confidence level, see Sec. 3.2), the corresponding uncertainties in H2 column densities 
at J = 2 and 3 can be relatively large. Fortunately, the physical state in the H2 absorber is 
not sensitive to Nj=2, 3 (see Sec. 4.1 for more details). Therefore we conclude that our final 
results are not critically dependent on the uncertainties in continuum fitting. 



3.4. Excitation temperature 

Once the H2 column densities at different J states are determined, we can calculate the 
excitation temperature for each J > 1 state by the Boltzmann law 

Nj gj r ej-eo, , . 

— = — exp TPf^l 3 



where ej—Eq is the excitation energy of level J relative to the ground state, k is the Boltzmann 
constant, the degeneracy is gj = (2J+1) x (2J+1), and / is the nuclear spin (0 for even J and 
1 for odd J). Fig. 6 shows the plot of log^ for each J state as a function of its excitation 
potential, with 1 a errors. The populations of the five detected J states can be well fitted by 
a single temperature Boltzmann distribution, with xl = 0.20. The one-temperature model is 
shown by the solid line in Fig. 6, with a best-fit excitation temperature of Tex = (152±10)K. 
The upper limit put on Nj=^ is also consistent with this model. Increasing or decreasing the 
best- fit continuum by 15% results in a similar excitation model with the best- fit excitation 
temperature of Tex = (153 ± 10) K or Tex = (150 ± 10) K, respectively. The insensitivity of 
Tex to the continuum fitting partly results from the logarithmic dependence of Tex on the 
ratio of column densities. 

A unique temperature for rotational excitation of H2 is not characteristic of Galactic 
H2 clouds, and the populations of rotational states of H2 often require at least two different 
excitation temperatures (e.g. Spitzer & Cochran 1973). Tex for the low J states approximates 
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Fig. 4. — The curve of growth for H2 hnes identified in the spectrum of Q1331+170. Different 
symbols stand for transitions from different J states. The background dots bracket the full 
range of the curve of growth for H2 lines, taking into account variations in the atomic data. 
The J — 0, 1 lines are on the damped part of the curve of growth, the J = 2, 3 lines on the 
flat part, and the J — 4 lines on the linear part. 



-23- 



19.4 - 



19.2 - 

o 



19.0 



O 



O 



o 



() 



C) 



C) 



J-0 



19.6- 
^ 19.4- 



O 



1 9.2 - 
18.8 ^ 



E 18.4 
^ 18.0 
9 17.6 



19.0 ^ 
^ 18.5- 

n 18.0 - 



17.5- 



1 5.t 



o 15.2 



o 



14.8 



14.4 



C) 



C) 



o 



C) 



o 



C) 



C) 



o 



o 



o 



o 



C) o 



o 



J = 2 



o 



o 



C) 



o 



o o 



o 



J = 3 



o 



o 



o 



C) 



o 



J = 4 



0.80 0.90 1.00 1.10 

Scaling of the best-fit continuum 



1 .20 



Fig. 5. — The effect of continuum fitting on the measured column densities of H2. 
abscissa shows the scahng factor apphed to the best-fit continuum shown in Fig. 1. 
best-fit continuum corresponds to a scahng factor of 1.0. 
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Fig. 6. — Distribution of molecular hydrogen at different J levels. The abscissa gives the 
relative energy for each J state relative to J = 0, expressed in units of Kelvin; the ordinate 
gives the column density divided by the degeneracy of the corresponding state, with la errors. 
The population of H2 at different J states can be fitted by a one-temperature Boltzmann 
distribution, with Tex ~ 152 K. 
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the kinetic temperature of the gas, while Tex for high J states results from UV pumping 
(Spitzer et al. 1974). For the damped Lya absorber of Q1331+170, a unique excitation 
temperature for all states from J = to 4 indicates that the rotational populations of H2 
arc probably dominated by collisional excitation in a gas with kinetic temperature around 
150 K. Moreover, the local UV radiation field does not significantly populate the high J 
states of H2 molecules. We quantify the UV radiation field below (see Sec. 4.1). 

4. Discussion 

4.1. UV radiation field 

A simple analysis on the rate equation of the J = 4 state can yield an estimate of 
the UV photo-absorption rate of H2. Collisions do not provide an important population 
source for J > 4 states (e.g. Browning et al. 2003), unless the gas phase temperature is 
comparable with or higher than T ^ > 1700 K. However, such a high temperature 

can be excluded for the damped Lya absorber toward Q1331-I-170, either by the distribution 
of H2 at different J states (sec Fig. 6), or by the observation of 21 cm absorption line (Wolfe 
& Davis 1979, Chengalur & Kanekar 2000). 

Previous work shows that the J = 4 state is mainly populated by direct formation 
into this level and pumping from J = 0, while the dominant depopulation mechanism is 
spontaneous decay (e.g., Dalgarno & Wright 1972, Jura 1975). Assuming a steady state, 
these effects can be combined into the following equation 

P4,oi?absn(H2, J = 0) + 0.19i?dustn(HI)n(H) = A42n(H2, J = 4), (4) 

where 7?abs is the photo-absorption rate, n(Hl) is the neutral hydrogen number density, 
n(H) pa n(HI) + 2n(H2) is the total hydrogen number density, -Rdust is the H2 formation 
rate on the surface of dust grains, p4,o = 0.26 is the UV pumping efficiency from J = to 
J — A (Jura 1975), and A42 = 2.75 x 10~^s~^ is the spontaneous transition probability from 
J = 4 to J = 2 (Wolniewicz et al. 1998). For simplicity, we assume 19% of all H2 formation 
results in the population at J = 4 (Jura 1975), although more recent work suggests that this 
fraction may vary for different types of dust grains (Takahashi & Uehara 2001). 

The equilibrium between the formation and photo-dissociation of H2 can be written as 

i?dustn(HI)n(H) ^ 0.11i?absn(H2), (5) 

where we assume that approximately 11% of photo-absorption leads to photo-dissociation 
(Jura 1974). Combining Eqn. 4 and 5, we get i?abs = (7.6±2.4) x 10"^^s"^ We note that the 
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photo-absorption rate derived this way is independent of the neutral hydrogen density, since 
the relevant term i?n(HI)n(H) cancels from Eqn. 4 and 5. Therefore R^hs is independent of 
the adopted H I column density (see Sec. 3.5). 

The photo- absorption rate, Rahs derived above is based on the H2 population at the 
J — A state. Since J — A hues are optically thin (see Sec. 3.3), the photo-absorption rate 
within the H2 cloud should be approximately the same as that outside the cloud. Therefore 
i?abs is also characteristic of the UV background intensity in the absorber's environment. 
With this in mind, we use i?abs derived above to estimate the ambient UV radiation field at 
1000 A from the expression, 

it:abs = 4.0 X 10-105,hield J^i^ S-\ (6) 

1000 A,0 

where J ° ft; 3.2 x 10~^° ergs s"^ cm~^ Hz~^ Sr~^ is the UV intensity at 1000 A in the 

1000 A,© ° 

solar vicinity (Hirashita & Ferrara 2005), and -Sghieid corrects for the the H2 self-shielding 
and/or dust extinction. Since the J = 4 hues are optically thin, the self-shielding of H2 can 
be reasonably ignored, and we estimate the importance of dust extinction following Hirashita 
& Ferrara (2005), i.e. 

-^shield ~ exp [-0.879ac( J ], (7) 

10^^ cm""' 

where k is the dust-to-gas ratio, defined as the ratio of the extinction optical depth to neutral 
hydrogen column density in units of 10^^ cm~^ (Pei et al. 1991). Vladilo (1998) showed that 
K can be estimated by 

« = 1.7 X 10[^"/"1(1 - 10^^/^"]). (8) 

Assuming that Znll and Fell are the dominant ionization states for these two species, we 
adopt [Fe/H] = -1.22 and [Fe/Zn] = -0.87 for the Q1331+170 DLA from Prochaska & 
Wolfe (1999). This gives k ^ 0.088, and therefore S'sMeid ~ 0.89, with our adopted H I 
column density (see Sec. 3.2). Inserting the above values into Eqn. 6, we get 

J ° 

loooA ^ (2.1 ±0.7) X 10-^ (9) 

1000 A,0 

i.e. the ambient UV radiation field of the H2 cloud associated with DLA 1331-1-170 is about 
three orders of magnitude weaker than that in the Solar vicinity. Assuming the spectrum of 
the UV radiation field follows Ji, oc we estimate the UV intensity at the Lyman limit 

as 

■^9121 ~ ^^-^ ^ ^"^^ ^ 

This value is consistent with the mean metagalactic UV background intensity of J -^^^ ~ 
'^■^^aio ^ 10~^^ ergs s~^ cm~^ Hz~^ Sr~^ at similar redshifts determined from proximity effect 



-27- 



in the Lya forest (Scott et al. 2002), implying that the ambient UV radiation of the H2 
absorber may originate externally from the cloud. If instead the UV field is stellar in origin, 
"^912 A lower than the metagalactic value, so our assumption of oc z/~°'^ gives 

a reasonable upper limit to the ionizing radiation field. Based on the above calculations, 
the existence of any significant internal source of UV emission is not allowed, suggesting 
negligible star formation activity associated with this H2 absorber. Finally, we emphasize 
that our estimate of J^^^ ^ is based on the assumption of QSO-only spectral shape for the 
UV background intensity. This is a reasonable assumption if star-forming galaxies do not 
make significant contribution to the metagalactic UV radiation field (e.g. Leithere et al. 
1995, Hurwitz et al. 1997, Deharveng et al. 2001, Giallongo et al. 2002, Fernandez-Soto et 
al. 2003). 

The ambient UV intensity for DLA 1331+170 is much smaller than that in other DLAs, 

which is usually close to the typical Galactic value (e.g. Ge & Bechtold 1997, Levshakov et al. 

2002, Ledoux et al. 2002). Hirashita & Ferrara (2005) estimated that the probable range of 

J 

radiation field is 0.5 < log ( j ) < 1.5, for a sample of DLAs which have been previously 

1000 A, 

searched for H2 absorption (Ledoux et al. 2003). The associated surface star formation rate 
(SFR) is 5 X 10-=^ Mq yr-i kpc'^ < Ssfr < 5 x lO'^ Mq yr-^ kpc'^ (Hirashita & Ferrara 
2005, see also Wolfe et al. 2003a, 2003b). Also, we note that Reimers et al. (2003) reported 
the detection of H2 associated with a sub-DLA system, with local UV radiation field inferred 

to be ~ 300 times stronger than the mean Galactic value. Although most DLAs may be 
associated with considerable on-going star formation, the DLA toward Q1331+170 at least 
represents a special system, which is in a stage cither before the onset of considerable star 
formation, or after star forming activities have ceased. 

For comparison, wc note that the C II* 1336 A absorption feature is useful for probing 
the ambient UV intensity near 1500 A (Wolfe et al. 2003a, 2003b). Thus it is interesting to 
compare the UV radiation field inferred by C II* with that determined from H2. Using the 
C II* absorption line as a diagnostic of the UV radiation field follows the idea that [CII] 158fim 
emission is the dominant coolant in DLAs, and is balanced by the grain photoelectric heating 
in a steady state. Since the grain photoelectric heating rate is proportional to the ambient 
UV intensity, a measurement of the strength of C II* 1336A absorption provides an estimate 
of the UV radiation field at 1500 A. Following Wolfe et al. (2003a), we express the heating 
rate per H atom as 

r = iO"''«^e(T7Pi^ -^''°°^-2u -ic -1 ) e^g« (11) 

10 ergs s ^ cm ^ Hz ^ Sr ^ 
where k, — 0.088 is the dust-to-gas ratio determined from Eqn. 8, e 0.049 is the heating 
efficiency (Bakes & Tielens 1994, Wolfire et al. 1995), and J-^^^^^^ is the UV intensity at 
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1500 A. For the DLA towards Q1331+170, the C II* 1336A absorption hne strength imphes 
a coohng rate of Ic ~ 2.24 x 10"^'' ergs H^^ (Wolfe et al. 2003a). Equating F and Ic imphes 
J " K, 5.4 X 10~^° ergs cm~^ s~^ Hz~^ Sr~^, not too far from the mean Galactic value and 

1500 A ^ ' 

much higher than the UV radiation field inferred from H2 absorption. This discrepancy is 
not expected, since the large difference in the UV intensity determined from H2 and C II* 
is difficult to be interpreted by any local sources of absorption. Note that H2 and C II* 
absorption probe UV intensity near ~ 1000 A and ~ 1500 A, respectively, which are both 
above the Lyman hmit. This indicates that photoionization of H I is not relevant. We also 
mentioned above that J ° determined from H2 absorption is characteristic of the UV 

1000 A z J- 

radiation field outside the H2 cloud, therefore molecular hydrogen absorption in the Lyman 
and Werner bands docs not account for the discrepancy either. However, the very low 
ambient UV intensity determined for the H2 absorber is reasonable, otherwise the residual 
excitation temperature inferred from C I/C I* absorption would be far below the cosmic 
microwave background (CMB) temperature predicted by the standard cosmology (see Sec. 
4.7). 



4.2. The structure of the H2 absorber 

To study the structure of the II2 cloud associated with the DLA of Q1331 + 170, we use 
CLOUDY (version c9400) to construct a grid of models with different total hydrogen number 
density n(H) k, n(HI) + 2n(H2), and examine which model provides the best match with the 
observed H2 column density. We adopt a metallicity of [Zn/H] = —1.22 (Prochaska & Wolfe 
1999) for all the elements in the absorber. To take into account differential dust depletion, 
we determine the dust-to-gas ratio k according to Eqn. 8. The ambient UV intensity is 
characterized by the metagalactic UV radiation field, assuming a QSO-like spectral shape, 
as implied by results in Sec. 4.1. For simplicity, we assume plane parallel geometry with 
a uniform total hydrogen number density. The cloud is illuminated from both sides in our 
model, and we run the calculation until the neutral hydrogen column density reaches the 
observed value of A^(HI) = 1.5 x 10^^ cm~^. 

The model with a total hydrogen number density, n(H) ^ n(III) + 2n(H2) ~ 0.22, gives 
the best match to the observed H2 fraction for the z = 1.7765 DLA system. We show in Fig. 
7.1 the thermal structure of the cloud, with the electron temperature determined from the 
balance between cooling and heating. In Fig. 7.2, we give the structure of the cloud traced 
by H I, H II, H2 and e. We find that the neutral hydrogen number density remains almost 
constant from the boundary deep into the central region. The molecular hydrogen fraction 
is very low in an outer shell with a thickness of ~ 46 pc, and increases sharply from < 10~^ 
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Fig. 7. — Model calculations of the thermal and ionization structure of the H2 cloud (see 
Sec. 4.2). The electron temperature (upper panel), the number densities of H I, H II, H2, 
electrons (middle panel), and the fraction of different species of carbon (lower panel) are 
each plotted as a function of the depth into the cloud. Note that only half of the cloud is 
shown in the figure (~ 1.2 kpc). 
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in the shell to ~ 10~^ in the center. We note that the outer shell is very thin, compared with 
the depth of the center region (~ 1.2 kpc). The ^^^^ ratio decreases monotonically with the 
depth into the cloud, and the metagalactic UV radiation field does produce some fraction 
of ionization in the outer shell of the cloud (^py ~ 20%). Fig. 7.1 shows that the electron 
temperature in the center H2 core is much lower than that in the outer shell, and has a mean 
value of Te ~ 140 K, in broad agreement with the excitation temperature determined in Sec. 
3.4. 

Fig. 7.3 shows the ionization structure of carbon for our model. The dominant ionization 
state is C 1 and C 11. Although the ambient UV radiation results in some fraction of C III 
before it's shielded by molecular hydrogen, the model calculation shows that the column 
density of highly ionized carbon species (e.g. C IV) is completely negligible. 

Compared with observations, these results suggest that the Q1331 + 170 absorber has 
multiple velocity components with different carbon ionization. First, the observed value of 
Ncn/Nci ^ 2647 for Q1331+170 (Liszt 2002) is much larger than the model prediction of 
~ 3.5 for the H2 cloud. Second, there is strong CIV absorption associated with this damped 
LycK system, with A^(CIV) > 10^^ cm^^ (Prochaska & Wolfe 1999). Therefore the bulk of the 
observed C II and C IV are associated with other components, either in a warm/hot phase, 
or embedded in a UV radiation field much stronger than that derived for the H2 absorber 
(see also Sec. 4.1). 

Dessauges-Zavadsky et al. (2004) reported multiple velocity components associated with 
the DLA of Q1331+170, and they did find non-neghgible differences between the component 
structure of low-ion transitions and that of intermediate-ion transitions. We note that the 
H2 cloud identified in our spectrum is different from the five intermediate-ion velocity com- 
ponents by 56 km s""*^, 28 km s^-*^, 15 km s^"*^, —15 km s^-*^, and —34 km s~^ in the rest frame, 
respectively. On the other hand, it agrees with one low-ion component with detected C 1 
absorption, with a velocity difference of only 1.6 kms^^ in the rest frame. Songaila et al. 
(1994) reported two C I absorption features associated with the same system, which are 
different from the H2 cloud by 1.4 km s~^ and 19 km s~^ in the rest frame. The 1 a uncer- 
tainty is 1.2 kms~^, caused by wavelength calibration and line centroiding (see Sec. 3.2), 
and therefore we conclude that the molecular hydrogen absorption identified in our STIS 
spectrum is consistent with one low-ion component reported by Dessauges-Zavadsky et al. 
(2004), as well as one C 1 component reported by Songaila et al. (1994). This is as expected 
since H2 tends to be co-spatial with neutral carbon and/or other low ion species. 
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4.3. Mciss and dynamical state of the H2 cloud 

Based on the model shown in Fig. 7.2, we estimate the mass of the entire cloud as 
M ^ 1.36 X |7ri?3[n(HI)mH + n{B.2)mYi.,] ^ 6.5 x 10^ Mq, where we take R ^ 1.2 kpc, 
n(HI) 0.17 cm^'^, n(H2) ~ 0.05 cm^^, and the factor of 1.36 corrects for the abundance of 
hehum (Dickman 1978). We refer to this mass as Mlte, since our model is constructed based 
on the assumption of local thermodynamic equilibrium. Here we have ignored the small mass 
fraction contributed by H II. We also assume spherical geometry for the absorbing region. 

The masses of giant molecular clouds (GMCs) have been measured for the Milky Way 
(e.g. Solomon et al. 1987, Digel et al. 1996, Heyer et al. 2001), and nearby galaxies 
such as LMC (Mizuno et al. 2001a), SMC (Mizuno et al. 2001b), and M33 (Engargiola et 
al. 2003). These masses are either determined based on the assumption of self-gravitational 
equilibrium, i.e. virial mass, Mvir, or determined by adopting a constant CO-to-H2 conversion 
factor and then integrating the H2 column density over the projected area of the cloud, i.e. 
CO luminosity mass, Mqq. In all cases, the mass spectrum of GMCs has been found to follow 
approximately a power law, with the power index ranging from ~ —1.8 for the Milky Way 
(Heyer et al. 2001) to ~ —2.6 for M33 (Engargiola et al. 2003). The mass of the H2 cloud 
associated with the DLA towards Q1331-I-170 is much larger than the values determined 
for most GMCs in either the Milky Way or nearby galaxies. For example, it is almost two 
orders of magnitude greater than the value of 7 x 10^ M0 for the most massive CMC in M33 
(Engargiola et al. 2003). Similarly, if we adopt the CMC mass spectrum of ^ oc M~^'^ 
in the Milky Way, and take the limiting CO luminosity of 138 K km s^^ pc^, corresponding 
to 566 Mq (Heyer et al. 2001), we find that the fraction of GMCs with masses greater 
than 6.5 x 10^ Mq (the value derived for the H2 cloud associated with DLA 1331+170) is 
only 9 X 10"^ in the Milky Way This may simply result from our assumption of spherical 
symmetry, or from the fact that the spectral resolution of our STIS spectrum is not high 
enough to resolve individual molecular clouds in velocity space along the hne-of-sight. Note 
that the typical velocity dispersions measured for Galactic GMCs peak at ~ 1 km s~^ 
(Heyer et al. 2001), much smaller than the spectral resolution of ~ 10 kms~^ for our 
spectrum and the 6-parameter of ~ 20 km s~^ for the H I gas in DLA 1331+170. The total 
H2 mass inside a 2 kpc galactocentric radius has been determined for a sample of 17 nearby 
galaxies (Paglione et al. 2001), based on the observed CO emission and a standard CO-to-H2 
conversion factor of 1.6 x 10^° cm~^ km~^ s. These masses range from 1.6 x 10^ M0 to 
2.8 X 10^ Mq, much larger than our value of 6.5 x lO'^M©. This is consistent with the general 
deficiency of molecular hydrogen in damped Lya systems, especially when considering that 
the DLA towards Q1331+170 presents the largest H2 fraction in all high-z DLAs. 

It is interesting to compare the mass derived above, i.e. Mlte, with the mass predicted 
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by self-gravitational equilibrium, i.e. Myir. We estimate the virial mass by 

Mvir ^ = 9.6 X 10^ Mq. (12) 

Here is the one-dimensional velocity dispersion which we approximate as the 6-parameter 
determined in Sec. 3.2, R 1.2 kpc is the size of the cloud (see Sec. 4.2), G is the 
gravitational constant, and X is a dimensionless factor of order unity that depends on the 
geometry as well as the density profile of the cloud. We adopt K ^ 8.7 from Solomon et 
al. (1987). Correspondingly, the gravitational parameter, — ~ 14, indicating that 
the H2 cloud can not be bound by self-gravity. This is contrary to the situation for most 
Galactic GMCs, of which the self-gravitational equilibrium state has long been established 
(e.g. Solomon et al. 1987). More recently, Heyer et al. (2001) found that the dynamical 
state of Galactic GMCs varies with the cloud mass, in that Mqq > 10^ Mq clouds arc bound 
by self-gravity. One possible interpretation to the large aQ value is that the H2 cloud towards 
Q1331-I-170 is in hydrostatic equilibrium, i.e. bound by the pressure of an external medium. 
Here we emphasize that the above result is subject to uncertainties due to the unknown 
geometry and density profile of the cloud. 



4.4. Moleculctr hydrogen fraction 

In Sec. 3.2, we derived a molecular hydrogen fraction of /hj = (5.6 ± 0.7)%, for the 
z — 1.7765 DLA toward Q1331-I-170, which is the largest value reported so far in any 
redshifted damped Lya absorber. The large value of /hj is consistent with the relatively 
high dust abundance for this system. In Sec. 4.2, we derived the dust-to-gas ratio to be 
k ~ 0.088. Ledoux et al. (2003) showed that there exists a threshold at A; ~ 0.03 above 
which the molecular hydrogen fraction tends to be large (/hj > 10~^). Moreover, the model 
described in Sec. 4.2 suggests a low gas phase temperature, which is also supported by the 
distribution of H2 at different J states. The low kinetic temperature is representative of a 
cold neutral medium (CNM) phase, in which efficient formation of H2 molecules onto dust 
grains is possible. Third, the ambient UV radiation field of the H2 cloud is extremely low, 
implying that destruction of H2 molecules by photo-dissociation is ineffective. 

It is interesting to compare the damped Lya absorber toward Q 133 1-1-1 70 with other 
high-redshift damped systems. In Table 4, we compiled data for all the 42 DLAs which 
have been searched for molecular hydrogen absorption, including Q1331-I-170. In the upper 
panel of Fig. 8, we show the distribution of damped Lya absorbers with respect to neutral 
hydrogen column density and molecular hydrogen fraction. The figure clearly shows the bi- 
model distribution of DLAs, i.e. for a certain value of A^hi, the molecular hydrogen fraction 
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Fig. 8. — Relation between molecular hydrogen fraction /hj and neutral hydrogen column 
density (top) and dust abundance (bottom) for damped Lya absorbers. The upper panel 
shows the bi-model distribution of high-redshift DLAs, which can be understood as a result of 
self-shielding of molecular hydrogen. The lower panel shows a strong correlation between the 
molecular hydrogen fraction and the dust abundance, characterized by the ratio of undepleted 
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of the absorber can be either as large as ~ 10~^ — 10~^, or smaller than the detection limit 
(~ 10~^ — 10~^). This can be interpreted by the self-shielding effect of molecular hydrogen: 
once molecular hydrogen begins to form, it shields itself from subsequent photodissociation 
by UV photons, and thus the molecular hydrogen fraction increases significantly (sec also 
Ge & Bechtold 1999). The lower panel of Fig. 8 shows the relation between the molecular 
hydrogen fraction and the dust abundance in the same system. The dust abundance is 
characterized by the column density ratio of undepleted element to depleted element. For 
element which is not depleted, we take Zn as a representative, whereas when Zn is not 
available in the literature, we adopt S instead. The element which is depleted to the dust 
grain is represented by Fe, and when Fe abundance has not been measured, we take Cr 
instead. The lower panel of Fig. 8 shows a strong correlation between the dust abundance 
and the molecular hydrogen fraction, similar to previous results (e.g. Ledoux et al. 2003). 
Kendall's r test shows the correlation is present at > 99.9% confidence level. In Fig. 8, the 
damped Lya absorber toward Q1331H-170 is indicated by a solid circle. 



4.5. Dust extinction 

The reddening of QSOs with intervening DLAs compared with non-DLA QSOs has been 
detected and quoted as evidence for the presence of dust in DLAs (e.g. Fall & Pei 1989, Pei 
et al. 1991). More recently. Murphy & Liske (2004) showed that there is no evidence for the 
dust reddening of QSOs with DLAs, based on a much larger sample selected from the Sloan 
Digital Sky Survey Data Release 2 (SDSS DR2) catalog. However, since the SDSS main 
QSO sample is largely color-selected (Richards et al. 2002), the non-detection by Murphy 
& Liske (2004) may simply be a selection effect against heavily reddened objects. In this 
section, we estimate the dust reddening of Q1331-I-170, due to the damped Lya absorber at 
z ~ 1.7765, and we also compare the reddening of QSOs for which the intervening DLAs 
have different molecular hydrogen fractions. 

Similar to Pei (1992), we calculate the color excess, -Bb-v, of Q1331-I-170 by 

where k ^ 0.088 for DLA 1331 + 170 is the dust-to-gas ratio determined from Eqn. 8, and 
R\ is the ratio of total-to-selective extinction, which depends on the exact extinction curve 
adopted in our calculation. If we assume the SMC extinction curve, then — 2.93 (Pei 
1992) and E^-w ^ 0.037 ± 0.005. We also calculate Eb-v by adopting the mean N^ji/E^-y 
values of 36.1 ± 3.3 based on the new determination of the SMC extinction curve (for the 
SMC bar sample) by Gordon et al. (2003), where A^hi is in units of 10^^ cm~^. This gives 
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the same value of -Bb-v within 1 a errors. There are two reasons for adopting the SMC 
extinction curve. First, the DLAs are beheved to be in an early stage of chemical evolution, 
which is well characterized by the SMC extinction curve (e.g. Pei 1992). Second, DLAs show 
no evidence for the 2175 A dust feature observed in either Galactic or LMC extinction curve 
(e.g. Fall & Pei 1989, but see Malhotra 1997 which reported the detection of this feature 
based on a sample of Mg II absorbers). 

In column 7 of Table 4, we give i?B-v for each QSO of which the intervening DLA has 
previous measurements of both depleted and undepleted metal absorption. These values are 
calculated similar to the Q1331+170 case. However, wc use a different normalization when 
estimating the dust-to-gas ratio by Eqn. 8, depending on the types of depleted/undepleted 
elements adopted in the calculation (see Table 1 of Vladilo 1998). For several QSOs with 
more than one intervening DLAs (Q0013-004, Q0405-443, Q0841+129 and Q2059-360), we 
only list in Table 4 £^b-v contributed by individual absorbers, which does not represent the 
total extinction towards the background QSO. Fig. 9 presents the distribution of i?B-v for 
DLAs with firm II2 detections (solid line) as well as the distribution for DLAs for which only 
upper limits of /hj have been put (dashed line). The extinction for Q1331+170 is marked 
with an arrow. The two distributions in Fig. 9 are different at 99.8% confidence level, based 
on Kolmogorov-Smirnov test. Specifically, the sub-sample of DLAs with firm H2 detections 
shows more extinction compared with other DLAs. The mean color excess for absorbers 
with firm H2 detections is (Eb-v) — 0.038 ± 0.001, significantly greater than the value of 
(£'b_v) = 0.014 ± 0.001 for cases in which only upper limits of have been reported. 
Whereas only 12.5% of DLAs with firm H2 detections have -Bb-v smaller than 0.01, this 
fraction is as large as 75% for other DLAs. We have also examined the relation between 
£^b-v and /hj, but no significant correlation has been found. 

The intrinsic QSO continuum in the UV/optical band is reddened by the dust extinction 
due to an intervening DLA. Assuming the QSO continuum follows F^, oc or equivalently 
F\ oc A"~^, wc get a 4.1 by fitting a power law to the dereddened continuum (between 
2600 A and 3000 A) of the QSO. Q1331 + 170, while a = 4.2 for the observed, reddened 
continuum adopted in Sec. 2. For dcrcddcning of the Galactic extinction, we use the 
extinction map of Schlegel et al. (1998). For reddening of the intervening damped Lya 
absorber z 1.7765, we adopt the FM parametrization (Fitzpatrick & Massa 1990) 
determined for the SMC bar sample (Gordon et al. 2003). The observed spectrum has 
been shifted to the rest frame of the absorber before dereddening. The value of 4.1 for the 
dereddened spectral slope indicates an extremely red intrinsic spectrum below the rest frame 
wavelength of ~ 970 A for Q1331+170. 
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Fig. 9. — The distribution of -Eb-v for damped Lya absorbers which have been searched for 
molecular hydrogen absorption. The solid line represents the histogram for the eight DLAs 
with firm H2 detections, whereas the dashed line represents the histogram for other DLAs 
for which only upper limits of /hj have been put. The two distributions are clearly different, 
in the sense that DLAs with H2 detections show significantly more extinction. 
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4.6. CO-to-H2 ratio 

According to our model calculations in Sec. 4.2, the predicted CO column density is 
6.7 X 10^°cm~^. This is consistent with the observed 2a upper limit of A^'co < 1-1 x lO^^cm"^ 
(Levshakov et al. 1988), determined from UV absorption lines of CO [A^Ii X^E^). 
Adopting the total molecular hydrogen column density in Table 3, we obtained the observed 
CO-to-Hs column density ratio of ^ < 2.5 x 10"^. 

The upper limit of the CO-to-H2 column density ratio for DLA 1331+170 is similar 
to the typical value of CO-to-H2 measured from UV absorption of CO in Galactic diffuse 
clouds, where the effects of CO self-shielding are not important (e.g. Crenny & Federman 
2004). For dense molecular clouds in the Milky Way, the CO/H2 conversion factor is given 

in terms of the quantity where Ico is related to the antenna temperature Ta directly 

H2 

measurable from miUimeter CO emission hnes (e.g. Strong et al. 1988). usually does 
not provide information on the corresponding column density ratio, since the profiles of 
the optically thick CO lines give the velocity width of molecular clouds rather than the CO 
column densities. The ratio of CO to H2 column densities in typical dense molecular clouds in 
the Milky Way is ~ 10^^, three orders of magnitude higher than the upper limit determined 
for the DLA toward Q1331+170 (e.g. Frerking et al. 1992, Lacy et al. 1994). This could be 
associated with the relatively low number density of the H2 cloud (see Sec. 4.2), since the 
CO abundance is controlled by coUisional reactions involved with various species, including 
O, C, H2, CH and CH2 (e.g. Tielens & HoUenbach 1985). 



4.7. Measurement of CMB Temperature at 2; = 1.77654 

Songaila et al. (1994) reported the detection of C I/C I* absorption features associated 
with a velocity component at 2; = 1.77654. As described in Sec. 4.2, the redshift of this 
component agrees with the redshift of the H2 absorber identified in our spectrum at 1.2 cr. 
Therefore, the C I and H2 probably arises in the same gas cloud. Theoretically, C I and 
H2 tend to be co-spatial as well (Ge et al. 1997, 2001). With the physical parameters 
determined in Sec. 4.1 and 4.2, we can estimate the local contribution to the C I excitation 
in this system and then put constraints on the CMB temperature at the absorber's redshift. 
The equilibrium between the population and de-population of the C I ^Pq and ^Pi fine 
structure can be written as 

Nq{BoJ, + Toi + X! ^01%) = ^1 (^10 + Bwlu + Tio + X! ^io%) , (14) 
3 3 
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where i?oi and Rio are coUisional excitation and de-excitation rates, with representing 
different collision partners (H, He, e, p, or H2), Fqi and Fio are UV pumping rates, Aiq = 
7.93 X 10~^ s~^ is the probability of the C I fine structure transition ^Pi —^^ Pq (Bahcall & 
Wolf 1968). B^ily and BiqI^ represent excitation or de-excitation due to the absorption of 
ambient microwave photons, which can be expressed as 

BoJ. = SBioI. = f^^llZ] , s-\ (15) 
exp (23.6 K/Texj — 1 

The coUisional rates, i?oi R{q are taken from Launay & Roueff (1977), Johnson et al. 
(1987), Roueff & Le Bourlot (1990), Staemmler & Flower (1991), and Schroder et al. (1991), 
and are related through the Milne relation, 

i?io = ^< exp (23.6 K/Tk). (16) 

We take number densities of different species from the model described in Sec. 4.2. 
The various collision terms arc calculated to be i?™n(HI) = 8.6 x 10~^° s~^, R^^niYie) = 

3.3 X 10-13 g-i^ R^^n{}i2) = 4.2 x 10"" s-\ RlMp) = 2.3 x lO'^^ g-i^ ^^^^ RoMe) = 

1.4 X 10^1^ s^i, assuming a kinetic temperature of 152 K. We determine the UV pumping 
rates for C I by scaling from the mean Galactic values (e.g. Ge et al. 1997), and get 
Foi = 8.7 X 10~i3 s~i and Fio = 2.8 x 10"^^ s^^ Inserting these values into Eqn. 14 and 15, 
and taking = 0.125 ± 0.042 from Songaila et al. (1994), we get = (7.2 ± 0.8) K 
at 2; = 1.77654. Note that the uncertainty in the calculated Tex is associated with the errors 
in the measured C I/C I* column densities. The excitation temperature, Tex, calculated 
above is consistent with the predicted CMB temperature of 2.725 x {1 + z) = 7.566 K at 
z — 1.77654, within 1 a. 



5. Conclusions 

The main results of this paper are summarized as follows: 

1. Using the NUV spectrum with spectral resolution of ~ 10 km s"^ obtained with the 
E230M grating of HST/STIS, we detected strong molecular hydrogen absorption associated 
with the z — 1.7765 damped Lya system toward Q1331-I-170. The total H2 column density 
is iVH2 = (4.45 ± 0.36) x 10^^ cm~^, determined from simultaneous Voigt profile fitting of 26 
H2 lines. 

27V 

2. The molecular hydrogen fraction was determined to be /hj = NTii+2Nii ~ (5.6±0.7)%, 
which is the largest value reported so far in any redshifted damped Lya system. This is a 
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combined effect of a relatively high dust-to- gas ratio, a low gas temperature, and an extremely 
low ambient UV radiation field. 

3. We detect rotationally excited transitions of H2, with J — O — A. The relative column 
densities of the J — 0, 1, 2, 3, 4 states can be fit with a single temperature of Tex = 152± lOK. 
This suggests that the derived excitation temperature represents the kinetic temperature of 
the gas, and that photo-excitation is negligible. 

4. The photo- absorption rate was estimated to be i?abs — (7.6 ± 2.4) x 10~^^ s~^, 
corresponding to an ambient UV radiation field of J^^^^ ^ 7x 10"^^ ergs cm~^ Hz""^ Sr~^. 
This is comparable with the metagalactic UV background intensity at the same redshift, and 
implies an extremely low star-formation rate in the DLA environment. 

5. We constructed a simple model to describe the structure of the H2 absorber, with 
a best-fit total hydrogen number density of ~ 0.2 cm~^ and an electron temperature of 
Tc ~ 140K. The total mass of the model cloud is ~ 6.5 x lO^M©, greater than masses of most 
GMCs in the Milky Way and nearby galaxies. We also find that this mass is considerably 
smaller than the virial mass of the cloud, indicating that the H2 cloud is not in a state of 
self-gravitational equilibrium. Our model is consistent with the population of H2 at different 
J states, as well as the observational upper limit put on the CO column density. The 
observed CO-to-H2 column density ratio is < 2.5 x 10~^, characteristic of the typical 
value measured for diffuse molecular clouds in the Galactic ISM. 

6. Adopting the SMC-likc extinction curve, we calculate the extinction of Q1331+170 
to be -Eb-v = 0.037 ± 0.005, due to the intervening damped Lya absorber at 2; ~ 1.7765. 
We also find that the extinction of QSOs by foreground DLAs with firm H2 detections is 
considerably greater than those for which only upper limits of /hj have been put. 

7. The redshift of the H2 absorber identified in our spectrum is consistent with a C I/CI* 
absorber reported by Songaila et al. (1994). We calculated the local contribution to C I 
excitation, including UV pumping and collisions. The residual excitation temperature was 
determined to be Tex = (7.2 ± 0.8) K, consistent with the predicted Tcmb of 7.566 K at 
z = 1.77654. 
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this work. We also thank R.F. Carswell and G. Ferland for making the Voigt profile fitting 
program VPFIT and the photoionization code CLOUDY publicly available. This work is 
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Space Telescope Science Institute, which is operated by the Association of Universities for 
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Table 3. H2 column densities", 6-parameter and redshift 



Parameter 


Measured value 


1 a error 


z 


1.776553 


0.000003 


b (km s~^) 


13.85 


0.52 


\og{Nj=o) 


19.144 


0.066 


log(iVj=i) 


19.467 


0.043 


log(A^j=2) 


17.901 


0.187 


l0g(A^j=3) 


17.674 


0.184 


l0g(A^j=4) 


15.101 


0.125 


l0g(A',„|.al) 


19.648 


0.035 



^Column densities are in units of cm ^. 
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Table 4: H2 absorption in redsliifted DLA and sub-DLA systems 



Object 

QOOOO - 263 3.390 

QOOlO-002 2.025 

Q0013 - 004 1.973 

Q0013 -004 1.968 

Q0058 - 292 2.671 

QOlOO+ 130 2.309 

Q0102 - 190 2.370 

Q0112 - 306 2.418 

Q0112 - 306 2.702 

Q0112 + 029 2.423 

Q0135 - 273 2.800 

Q0347 - 383 3.025 

Q0405 - 443 2.550 

Q0405 - 443 2.595 

Q0405 - 443 2.621 

Q0454 + 039 0.860 

Q0458 - 020 2.040 

Q0515 -441 1.151 

Q0528 - 250 2.811 

Q0551 - 366 1.962 

Q0841 + 129 2.374 

Q0841 + 129 2.476 

Q0935 + 417 1.373 

Q1037 - 270 2.139 

QllOl - 264 1.839 

Q1117- 133 3.351 

Q1157 + 014 1.944 

Q1223 + 178 2.466 

Q1232 + 082 2.338 

Q1328 + 307 0.692 

Q1331 + 170 1.776 

Q1337+ 113 2.508 

Q1337+ 113 2.796 

Q1444 + 014 2.087 

Q1451 + 123 2.469 

Q1451 + 123 3.171 

Q1946 + 766 2.843 

Q2059 - 360 2.508 

Q2059 - 360 3.083 

Q2138 - 444 2.852 



Zabs logiV(HI) 



log /h2 



Pel b 
H 



0.01 



21.41 ±0.08 
20.80 ±0.10 
20.83 ±0.04 

< 19.43 
21.00 ±0.10 
21.40 ±0.05 
20.85 ±0.08 
20.37 ±0.08 
20.15 ±0.07 
20.70 ±0.10 
20.80 ±0.10 
20.56 ±0.05 
21.00 ±0.15 
20.90 ±0.10 
20.25 ±0.10 

20.69 ±0.06 
21.65 ±0.09 
19.88 ±0.05 
21.20 ±0.10 
20.50 ±0.08 
20.90 ±0.10 
20.65 ±0.10 
20.40 ±0.10 

19.70 ±0.10 
19.35 ±0.04 
20.85 ±0.10 
21.70 ±0.10 
21.40 ±0.10 
20.90 ±0.10 

21.28 
21.18 ±0.04 
19.95 ±0.05 
20.85 ±0.08 
20.07 ±0.07 
20.30 ±0.10 
19.90 ±0.20 
20.27 ±0.06 
20.14 ±0.07 
20.85 ±0.08 
20.82 ±0.05 



< -6.98 

< -6.09 
-2.81 ±0.64 

> -2.36 

< -6.69 

< -5.30 

< -5.60 

< -5.72 

< -5.50 

< -6.32 

< -6.40 
-5.71 ±0.10 

< -6.59 
-2.44 ±0.23 

< -6.15 

< -4.10 

< -5.70 
-2.64 ±0.30 
-2.58 ±0.22 
-2.78 ±0.90 



< 
< 
< 
< 
< 
< 
< 
< 
> 
< 



-5.93 
-6.05 
-3.70 
-5.40 
-4.75 
-6.37 
-6.69 
-6.80 
-3.41 
-5.30 
-1.25 ±0.05 

< -5.37 

< -6.54 
-1.48 ±0.38 

< -4.70 

< -5.80 

< -5.05 

< -4.93 

< -6.64 

< -6.87 



-2.07 ±0.10 
-1.20 ±0.12 
-0.75 ±0.06 

> -0.73 
-1.42 ±0.11 
-1.62 ±0.06 
-1.73 ±0.14 



-1.14 ±0.13 * 
-1.29 ±0.17 * 
-0.98 ±0.09 
-1.17±0.16 
-1.02 ±0.12 

-1.03 ±0.12 
-1.19 ±0.09 

-0.54 ±0.05 
-0.78 ±0.12 
-0.15 ±0.09 
-1.52 ±0.14 
-1.52 ±0.11 
-0.80 ±0.14 
-0.26 ±0.11 
-0.82 ±0.14 * 
-1.28 ±0.13 
-1.32 ±0.10 
-1.63 ±0.11 

-1.21 
-1.22 ±0.04 



-0.60 ±0.15 



-1.76 ±0.09 * 
-1.65 ±0.12 * 
-1.48 ±0.05 



-2.04 ±0.09 
-1.25 ±0.11 
-1.51 ±0.05 

> -2.33 
-1.76 ±0.10 
-1.93 ±0.05 
-1.89 ±0.13 
-2.50 ±0.09 
-0.89 ±0.10 
-1.35 ±0.11 
-1.54 ±0.15 
-1.72 ±0.06 
-1.52 ±0.15 
-1.33 ±0.11 
-2.15 ±0.10 
-1.02 ±0.08 
-1.77 ±0.10 

-1.55 ±0.10 ** 
-1.25 ±0.15 
-0.95 ±0.09 
-1.71 ±0.11 
-1.51 ±0.11 

-0.90 ±0.20 ** 
-0.65 ±0.10 
-1.36 ±0.05 
-1.55 ±0.13 
-1.73 ±0.10 
-1.70 ±0.10 
-1.73 ±0.13 

-1.66 
-2.09 ±0.04 
-2.05 ±0.07 
-2.02 ±0.09 
-1.58 ±0.09 
-2.41 ±0.10 
-2.09 ±0.24 
-2.53 ±0.06 
-2.12 ±0.05 
-1.83 ±0.18 
-1.68 ±0.05 



0.00 ±0.33 
0.20 ±0.67 
4.65 ±0.90 

0.97 ±0.54 
1.44 ±0.45 
0.19 ±0.24 



0.47 ±0.43 
0.48 ±0.49 
1.46 ±0.41 
1.75 ±1.40 

1.81 ± 1.15 

0.00 ±0.72 
9.95 ±3.57 
0.96 ±0.28 
8.14 ±4.14 

8.82 ±2.73 
0.40 ±0.42 
0.00 ±0.23 
0.40 ±1.92 
0.76 ±0.39 
0.08 ±0.04 
0.81 ±0.62 
6.86 ±3.19 
0.41 ±0.89 

3.55 
3.69 ±0.52 



1.24 ±0.52 



0.02 ±0.17 
0.17±0.17 
0.38 ±0.15 



